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Phase equilibria of aqueous solutions of formaldehyde and methanol: improved 
approach using UNIQUAC coupled to chemical equilibria 
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a Universite de Toulouse; INPT, UPS; Laboratoire de Genie Chimique; 4, Allee Emile Monso, F-31030 Toulouse, France 
b CNRS; Laboratoire de Genie Chimique; F-31030 Toulouse, France 


Abstract 

The study of the phase equilibria involving formaldehyde is still relevant because of its presence in new processes 
where biomass is the raw material. The coupling between physical phase equilibria and chemical reactions makes 
its thermodynamic description a challenging task. In this work, an improved approach using UNIQUAC coupled to 
chemical equilibria was developed and compared with experimental data from the literature. The first application was 
done for the phase equilibria of the formaldehyde-water system and distribution of oligomers in the liquid phase was 
computed. The second and the third applications respectively considered the phase equilibria of the formaldehyde- 
methanol system and the formaldehyde-water-methanol system. 

Keywords: phase equilibria, modeling, formaldehyde, water, methanol, UNIQUAC 


1. Introduction 

Formaldehyde has been used for long as intermediate chemical in several industrial processes (such as polymers 
production [1], adhesive synthesis [2], trioxane production [3], formaldehyde distillation [4,5]...). Because of its acute 
toxicity (corrosivity, carcinogenicity, mutagenicity and reprotoxicity), its use is declining but studies about systems 
involving this compound are still relevant either because they are useful for depollution processes [6, 7] or because 
emerging processes where lignocellulosic biomass is the raw material are likely to generate this compound. This is for 
instance the case for thermal processes like torrefaction [8, 9] where formaldehyde was shown to be present in large 
amount in the gaseous effluent alongside other compounds including water and methanol. Because of its specific 
nature, modeling the behavior of formaldehyde is an important step for understanding and controlling of thermal 
processes involving lignocellulosic biomasses in general. Indeed, formaldehyde is also present in pyrolysis oils [10] 
and should be taken into account for depollution processes. 

Formaldehyde is the smallest aldehyde molecule and is a gas at ambient conditions. It is highly soluble and reactive 
in water. Therefore formaldehyde is commonly handled in aqueous and/or methanolic solutions that stabilize it, the 
most common being known as formalin, an aqueous solution of formaldehyde and methanol containing between 37% 
and 41% of formaldehyde. Indeed, aqueous solutions of formaldehyde and methanol are not simple ternary systems 
because formaldehyde reacts with both methanol and water to form diverse polymers. Thus, formaldehyde is not 
stored or processed as a pure substance. Reactions with water generate methylene glycol, and poly(oxymethylene) 
glycols. Reactions with methanol form hemiformal and poly(oxymethylene) hemiformals. 

This leads to a great complexity for the description of phase equilibria of this system and surprisingly, very 
few works are present in the literature. The most comprehensive studies originate from the group of Maurer at the 
University of Kaiserslautern (Kaiserslautern, Germany) and constitute the reference for this domain. Data from this 
group were used in this work [11, 12, 13, 14, 15, 16, 17, 18]. The model established by this group was used as a basis 
for our improved approach that was assessed by comparison to experimental data and the original model. 
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2. Description of the reactive vapor-liquid equilibrium model 

As mentioned above, when modeling the thermodynamic behaviour of such systems, the main difficulty is to 
account for the coupling of chemical and physical equilibria of these reactive molecules. In this work, the physical 
and the chemical phenomena were implemented in a uncoupled way in the model so as to differenciate the effects 
of weak intermolecular interactions of the physical equilibria from the strong intermolecular interactions involved 
in the chemical reactions. For the vapor-liquid equilibrium an heterogeneous approach was adopted: the physical 
interactions between all species are taken into account through activity coefficients calculation in the liquid phase and 
through an equation of state for the gas phase. In this study, chemical reaction equilibrium constants are not considered 
as variables. The only variables to be estimated are the binary interaction parameters. The figure 1 illustrates the 
outline of this model. Note that the system is described at equilibrium and no kinetic data are introduced in the 
modeling. Thus, the thermodynamic problem includes both: 

• chemical equilibria of the methylene glycol, hemiformal, poly(oxymethylene) glycols, poly(oxymethylene) 
hemiformals formation. 

• physical phase equilibria of water, methanol, formaldehyde, methylene glycol and hemiformal. 



Figure 1: Scheme of the vapor-liquid phase and chemical equilibria for aqueous solution of formaldehyde and methanol 

This description of phase equilibrium in chemical reactive mixtures was fomerly proposed by Maurer [11] and 
applied to aqueous solutions of formaldehyde and methanol. The UNIFAC Original model was chosen to represent 
the physical phase equilibrium. The advantage of the UNIFAC Original method lies in its predictive aspect but an 
important limitation is the influence of so-called proximity effects which are not accounted for. Indeed, for mixtures 
containing small molecules, the environment has a strong effect on the phase equilibrium. In this work, chemical de¬ 
scription of the Maurer’s approach was not modified but thermodynamic approach was improved using the UNIQUAC 
model. The main interest of the UNIQUAC model is the use of available experimental data for binaries to provide 
a more realistic description of the mixture behavior and for instance does not suffer from the limitations of UNIFAC 
Original, such as proximity effect. Moreover, it is able to account for size effects. Note that our approach is intended 
to be extended to the description of more complex mixtures containing other small polar molecules as encountered in 
gaseous effluents of wood torrefaction processes. A representative model of the system would allow proposing and 
designing a reliable separation process. In this case, the sole use of the UNIFAC Original model would give imperfect 
prediction. 

2.1. Chemical reactions in aqueous and methanol formaldehyde solutions 

Chemical equilibria are included in the thermodynamic description and chemical equilibrium is assumed. Formalde¬ 
hyde is a very reactive component. In this model, the prominent reactions were assumed to be: 
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• formation of methylene glycol (MG): CH 2 0 + H 2 0 ^ HO{CH 2 G)H 

• formation of poly(oxymethylene) glycols (MG„): H0{CH 2 0) n -\H + H0(CH 2 0)H =; HO{CH 2 C)) n H + H 2 0 

• formation of hemiformal (HF): CH 2 0 + CH 3 OH — CH 3 0(CH 2 0)H 

• formation of poly(oxymethylene)hemiformals(HF„): CH 3 0(CH 2 0) n -\H+CH 3 0(CH 2 0)H ;= CH 3 0(CH 2 0) n H+ 
CH 3 OH 

The more concentrated the formaldehyde solution, the higher the degree of polymerization (see figure 2). Maurer’s 
works considered polymers up to degree 4. Nevertheless, in this work, to obtain mass balance accuracy better than 
5%, it was necessary to 

100 % 

90% 

= 80% 

1 70% 
t 60% 

| 50% 

2 40% 

J- 30% 

a 20% 

H 10% 

0 % 

Overall liquid molar fraction of formaldehyde x FA overal1 


onsider polymers up to degree 7. 



6% 13% 19% 20% 21% 28% 29% 35% 39% 45% 


Figure 2: isothermal phase equilibrium for water-formaldehyde system at 353K [17]:species distribution vs overall molar fraction of formaldehyde 
in the liquid phase 


Chemical reaction equilibrium constants are taken from the literature and follow a polynomial law: InK = a\ + 
a 2 /T. The coefficients are given in table 1. 


2.2. Description of the liquid phase 

The major improvement proposed in this work is a better description of the thermodynamic behavior of the liquid 
phase. The physical interactions between all species are taken into account by an empirical approach based on local 
composition: the universal quasi-chemical model (UNIQUAC). UNIQUAC equations [20] are given by: 


In ji = In yf + In yf 


r Z 

In yf = In — + - In 

1 x i 2 

lny? = 9i(l-lnj]" =1 



Xjlj 


OkTkj 


( 1 ) 

( 2 ) 

(3) 


Z"=i ' Em *m 
= r° + tJj(T - T „ f ); t^ = t “ + TjfT - T ref ) 


(4) 

(5) 
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Table 1: Chemical reaction equilibrium constants: In K = ai + a 2 /T. 


Reaction 

Phase 

Heat of reaction (kJ/mol) 

ai 

a 2 

References 

W + FA & MG 

Vapor 

-43.51 

-16.984 

5233.2 

[19] 

2 MG o MG 2 + W 

Liquid 

-0.234 

4.98.10 -3 

869.5 

[17, 18] 

MG + MG 2 MG 3 + W 
MG + MG 3 o MG a + W 
MG + MG a o MG 5 + W 
MG + MG 5 o MG 6 + W 
MG + MG 6 o MG 7 + W 

Liquid 

-0.234 

1.908.10 -2 

544.5 

[17, 18] 

ME + FA ^ HF 

Vapor 

-53.73 

-14.755 

5969.4 

[14] 

2 HF o HF 2 + ME 

HF + HF 2 o HF 3 + ME 
HF + HF 3 o HF a + ME 
HF + HFa o HF 5 + ME 
HF + HF 5 o HF 6 + ME 
HF + HF 6 o HFj + ME 

Liquid 

-7.00 

-0.4966 

-491.3 

[14] 


UNIQUAC is an empirical model that requires experimental data to identify binary interaction parameters. Therefore, 
it is not totally predictive as the UNIFAC model. This model has proved to give a good description of polar mixtures. 
Also, it makes it possible to handle complex multicomponent mixtures, as it is the case here, from the knowledge of 
identified binary interaction parameters of all involved binary sub-systems. In the case of the formaldehyde-water- 
methanol system, the model is likely to describe the behavior of both small molecules (FA, W, ME, HF, MG) and larger 
ones (MG2 - MG7, HF 2 - HF 7 ). Indeed, as we already mentioned, the UNIQUAC model was also chosen for its ability 
to account for the influence of the size and the form of the molecules in a mixture thanks to the combinatorial term of 
equation 1. The residual term of equation 1 characterizes energetic interactions between molecules and requires the 
identification of the binary interaction parameters. The temperature dependency of the binary interaction parameters 
Tjj and Tji is also accounted by equation 5. The identification of the binary interaction parameters is developed in 
the section 3. The oligomer distribution in the liquid phase is obtained from mass balances, assuming chemical 
equilibrium. Overall compositions are calculated from mass balances. They are given by equations 6 to 9. 


x71 r jf = s -(XfA + XmG + 7 ,._ 0 i-XMGi + XHF + / ,._ 0 i-XHF;) 

(6) 

xf^ emI1 = S.(x W + XmG + XL'-I 

(7) 

X°m™ U = sfx ME + X HF + X HFj ) 

(8) 

S = (1 + XmG + y, : _ n i-XMGi + X HF + i-XHF,)~' 

(9) 


These fractions are calculated for every Pxy diagram. 

2.3. Description of the vapor phase 

From a physical point of view and as the model is to be used in low pressure conditions, the vapor phase is assumed 
to behave as a mixture of ideal gases. It contains water, methanol, formaldehyde, methylene glycol and hemiformal. 
The chemical reaction equilibria for the formation of methylene glycol and hemiformal were accounted for through 
the chemical-reaction equilibrium constants described in the next section. 
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2.4. Equations of the reactive vapor-liquid equilibrium model 

The uncoupled approach adopted to model the reactive vapor-liquid equilibrium leads to equations 10 to 15. The 
equations for the saturation vapor pressures (see equation 11) of the components were taken from literature [18] and 
their coefficients are given in table 2. 

y i P = x i y i (T,x)P s i (J) (10) 


In P s i {T)/kPa=A 


Kmg„ = 


T/K + C 

_ yMG P ref 

yFAyw p 
XMG n Xw 7mG„7w 
xmgXmg „-i 7 mg7mg„_i 
= y H F Pref 
yFAyME p 
XHF n XME jHF n yME 
XHFXHF'jHFjHF n . x 


(ID 

( 12 ) 

(13) 

(14) 

(15) 


As explained above, the activity coefficients are calculated in this work with the UNIQUAC model (see equations 1 
to 5). 

Note that when the vapor-liquid equilibrium equations and the chemical reaction equilibrium equations in one 
phase are satisfied, the chemical-reaction equilibrium equations in the other phase are automatically satisfied. 


Table 2: Antoine coefficients for pure component vapor pressures lnP s = A + B/(T + C) [18] 


Component 

A 

B 

C 

Formaldehyde 

14.4625 

-2 204.13 

-30.0 

Water 

16.2886 

-3 816.44 

-46.13 

Methanol 

16.5725 

-3 626.55 

-34.29 

Methylene Glycol 

19.5527 

-6 189.19 

-9.15 

Hemiformal 

19.5736 

-5 646.71 

0.00 


toe 3. Determination of the UNIQUAC binary interaction parameters 

10? The not reactive water-methanol binary has been widely studied in the literature. Nevertheless, as the reported 
io8 binary interaction parameters may have been estimated with other values of the pure component properties than those 
loo used in this study, they were identified again in this work using experimental data recommended by the DECHEMA 
no in the pressure and temperature range of [267 mbar-1013 mbar] and [35°C-100°C] respectively [21, 22], 
m Other binaries involve chemical reactions. To avoid spreading the uncertainty of the value of the chemical equi- 
112 librium constant into the physical equilibrium parameters, the binary interaction parameter estimation has been per- 
n3 formed on non-reactive VLE data. In consequence, only vapor-liquid equilibrium equations are needed to estimate 
n4 the binary interaction coefficients (equations 1 to 5 and 10 to 11). As it is impossible to uncouple chemical and 
ns physical equilibria for these systems in practice, no experimental data are available for the non-reactive vapor-liquid 
ns equilibrium. So simulated vapor-liquid data were generated using the UNIFAC Original model. 
n7 The objective function is usually defined as the residual between experimental data and estimated vapor-liquid 
ns flash calculations. Indeed, it was not possible to use experimental data in our case because of the presence of the 

119 chemical reactions. Therefore, vapor-liquid data were generated using the UNIFAC Original model as developed by 

120 the group of Maurer [18]. In consequence, the objective function was defined as the residual between UNIQUAC and 

121 UNIFAC Original estimated pressures and mole fractions: 

Fobj,bubble + Fobfdew £ 

r obi = - ( 16 ) 

nbubbie + n dew 
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nbubb 

V~! ^bubble '-D'l.C/A 

r obj,bubble = / r _, - 


c - y bubble \) i (|yhfhh 


ubble i\ (\ pbubb 

.calc y\ r UNL 


VieW\n buU e\tt ble 

and 


(J’ubb, 

yi,»A 


ybubble 

ibbley 




ybubble ^. 


p bubble 

(pbubble . pbubble\ 

i,UNIF AC " r r i,calc H 


_ yim. (fyutJNIFAC yfTalc^i (^%UNIFAC %lj)i d^LW/FAC _ 

= 2 j,-=1 ft™ + + 


(17) 

(18) 

(19) 


'1ie\\-,n dew \tr = (/l ™ NIFAC 2 + /x7alc)i ; yg" = Cwac+C). , ^ = (^W + O* (20) 

UNIQUAC binary interaction parameters have to be identified for systems including water, methanol, formaldehyde, 
methylene glycol, hemiformal, poly(oxymethylene) glycols (MG2 - MG7), and poly(oxymethylene) hemiformals 
(HF 2 - HF 7 ). Without any assumption, binary interaction parameters of 136 binaries should have been identified. 
Some preliminary calculations, validated by using UNIFAC Original, showed that it was not worth differentiating 
the binary interaction parameters of the polymers of the methylene glycol and poly(oxymethylene) glycol with other 
components except for formaldehyde. This consists in neglecting the size effect in the second term of the residual part 
of the activity coefficient for these compounds. Finally, 33 set of parameters were identified by solving the following 
optimization problem: 

min F obj (21) 


The identification was done using the Matlab function fminsearch which is a multidimensional unconstrained nonlin¬ 
ear method of minimization (Nelder-Mead) and which was coupled with the Simulis Thermodynamics Toolbox in our 


137 4. Results and Discussions 

138 4.1. UNIQUAC binary interaction parameters 

139 The UNIQUAC binary interaction parameters were obtained by minimizing the objective function defined by 
mo equation 16. For all data points in this work, the relative mean error was calculated. The value of the objective 
in function indicated a good estimation of the UNIQUAC binary interaction parameters. UNIQUAC binary interaction 

142 parameters that were determined in this work are given in tables 3 to 6. 

143 Table 3 reports the water-methanol UNIQUAC estimated binary interaction parameters. 


Table 3: Water - Methanol UNIQUAC estimated binary interaction parameters 


Component 1 

Component 2 r? 

tJj Tj t rj/ References 

Methanol (ME) 

Water (W) 156 

-369 0.91 0.20 [21,22] 


144 Table 4 shows the UNIQUAC binary interaction parameters of the formaldehyde-water system and table 5 reports 

145 the UNIQUAC binary interaction parameters of the formaldehyde-methanol system. The table 6 shows the binary 
us interaction parameters to add to deal with the water-formaldehyde-methanol ternary system. 

147 4.2. Comparison with literature data 

ms To check the validity of this approach, the reactive vapor-liquid model was solved. The activity coefficients were 

149 calculated by using either the UNIQUAC model with binary interaction parameters identified in this work (chemU- 

150 NIQUAC model) or the UNIFAC Original model (chemUNIFAC model). The results were compared to experimen- 

151 tal data taken from the literature for the formaldehyde-water system [17], formaldehyde-methanol system [23], and 

152 formaldehyde-water-methanol system [18]. 
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Table 4: Formaldehyde-Water system UNIQUAC estimated binary interaction parameters 


Component 1 

Component 2 

r° 

Tjj 

T° 

it 

Formaldehyde (FA) 

Water (W) 

967 

-169 

1.94 

-0.23 

Formaldehyde (FA) 

Methylene Glycol (MG) 

2714 

190 

2.21 

-2.74 

Formaldehyde (FA) 

Methylene Glycol 2 (MG 2 ) 

-274 

506 

-0.05 

-1.72 

Formaldehyde (FA) 

Methylene Glycol 3 (MG3) 

-22 

153 

-0.28 

-1.75 

Formaldehyde (FA) 

Methylene Glycol 4 (MG4) 

-54 

149 

0.12 

-2.36 

Formaldehyde (FA) 

Methylene Glycol 5 (MG5) 

26 

75 

0.09 

-2.59 

Formaldehyde (FA) 

Methylene Glycol 6 (MG 6 ) 

32 

32 

0.24 

-2.82 

Formaldehyde (FA) 

Methylene Glycol 7 (MG7) 

35 

31 

0.32 

-3.08 

Water (W) 

Methylene Glycol (MG) 

431 

-399 

0.16 

-0.25 

Water (W) 

Poly(oxymethylene) glycols (MGi> 2 ) 

65 

-299 

0.18 

0.76 

Methylene Glycol (MG) 

Poly(oxymethylene) glycols (MG ; > 2 ) 

977 

-945 

-1.56 

2.18 

Poly(oxymethylene) glycols (MG,> 2 ) 

Poly(oxymethylene) glycols (MG j>2,j*v 

140 

-103 

-0.0067 

-0.034 


Table 5: Formaldehyde-Methanol system UNIQUAC estimated binary interaction parameters 


Component 1 

Component 2 

r° 


r° 

t t 

Formaldehyde (FA) 

Methanol (ME) 

577 

-251 

0.29 

-0.4 

Formaldehyde (FA) 

Hemiformal (HF) 

-54 

76 

-0.46 

-0.72 

Formaldehyde (FA) 

Hemiformal 2 (HF 2 ) 

-37 

58 

-0.28 

-1.36 

Formaldehyde (FA) 

Hemiformal 3 (HF 3 ) 

-25 

31 

-0.02 

-1.87 

Formaldehyde (FA) 

Hemiformal 4 (HF 4 ) 

-3 

11 

0.14 

-2.26 

Formaldehyde (FA) 

Hemiformal 5 (HF 5 ) 

-18 

18 

0.32 

-2.61 

Formaldehyde (FA) 

Hemiformal 6 (HF 6 ) 

-19 

17 

0.42 

-2.88 

Formaldehyde (FA) 

Hemiformal 7 (HF 7 ) 

-19 

16 

0.55 

-3.16 

Methanol (ME) 

Hemiformal (HF) 

26 

-38 

-0.85 

0.13 

Methanol (ME) 

Poly(oxymethylene) hemiformal (HF,> 2 ) 

12 

0.24 

-0.10 

-1.36 

Hemiformal (HF) 

Poly(oxymethylene) hemiformal (HF,> 2 ) 

-24 

23 

1.47 

-1.47 

Poly(oxymethylene) hemiformal (HF 1 > 2 ) 

Poly(oxymethylene) hemiformal (HF ; > 2 ^,-) 

37 

-26 

0.43 

-0.47 


Table 6: Formaldehyde-Water-Methanol system UNIQUAC estimated binary interaction parameters 

Component 1 

Component 2 

r° 

if. 

r° 

t t 

Hemiformal (HF) 

Water (W) 

825 

-16 

-2.97 

2.17 

Hemiformal (HF) 

Methylene Glycol (MG) 

-291 

-16 

1.35 

0.63 

Hemiformal (HF) 

Poly(oxymethylene) glycols (MG,-> 2 ) 

501 

-97 

2.58 

-1.09 

Poly(oxymethylene) hemiformal (HF,> 2 ) 

Water (W) 

457 

985 

-2.00 

0.41 

Poly(oxymethylene) hemiformal (HF,> 2 ) 

Methylene Glycol (MG) 

-21 

663 

0.46 

1.00 

Poly(oxymethylene) hemiformal (HF,> 2 ) 

Poly(oxymethylene) glycols (MG / > 2j;t! ) 

878 

-94 

0.00 

0.00 

Methanol (ME) 

Methylene Glycol (MG) 

15 

-130 

-0.39 

0.03 

Methanol (ME) 

Poly(oxymethylene) glycols (MGi> 2 ) 

206 

-333 

-0.84 

-0.55 


7 


Page 7 of 15 





ACCEPTED MANUSCRIPT 


153 For all systems, the average overall deviation in the gas-phase composition and the average overall deviation of 

154 the pressure between experimental data and our work were calculated as: 




(I yFA,exp - yFAjnodeldl 


( 22 ) 


1 \rine* p (\yME,exp-yME,model\)l 


yMEJ 

(I Pexp ~ Pmc 


(23) 

(24) 


z n z 

Apmodel = _L y"“ p 

n exp ^'=1 

with model = chemUNIQUAC or chemUNIFAC and V / e [ 1; n exp ]y fa.i = ; y ME l = ^ = 

(Pe X p+Pmodel)l 

Concerning the accuracy of the experimental data from the literature, formaldehyde concentrations were deter¬ 
mined with a relative error of less than 2%. The temperature and the pressure were measured with an accuracy of ± 
0.1K and ± 0.5 kPa respectively [17]. 


161 4.2.1. Water-Methanol System 

162 The average deviation of the gas-phase composition and average deviation of the temperature were calculated for 

163 the water-methanol binary and are reported in table 7. 


Table 7: Average deviation of the gas-phase composition and average deviation of the temperature for water-methanol binary using UNIQUAC. 


Compound 1 

Compound 2 

Type of diagram 

AT (%) 

Ayi (%) 

References 

Water 

Methanol 

P = 1013 mbar 

0.18 

1.90 

[21] 

Water 

Methanol 

P = 666 mbar 

0.43 

3.58 

[22] 

Water 

Methanol 

P = 466 mbar 

0.84 

5.46 

[22] 

Water 

Methanol 

P = 266 mbar 

0.57 

6.37 

[22] 


165 4.2.2. Formaldehyde-Water System 

166 Table 8 compares the average overall deviation of the gas-phase composition and average deviation of the pressure 
le? for the water-formaldehyde system between experimental data from the literature [17] and the two models (chemUNI- 

168 QUAC and chemUNIFAC). Furthermore, the average deviation of the polymer distribution in liquid-phase, between 

169 experimental data [17] and the two models (chemUNIQUAC and chemUNIFAC), were calculated using equation 25 

170 (see tables 9 and 10). 

—p = J_ y'w (\xM Gn , exp -x MGnJnodel \)i (25) 

" n exp x—U-t x M g„j 

171 with n the degree of polymerization and model = chemUNIQUAC or chemUNIFAC and V / e [1; n exp ] x M g„j = 



73 The values of the deviations indicate that the chemUNIQUAC model provides a fairly good description of the 

74 formaldehyde-water system. Figure 3 shows the isothermal phase diagram at 383K. Good results were obtained for 

75 the prediction of the overall composition (a) and the distribution of the polymers in the liquid phase (b). Our model 

76 (solid line) was compared with the original chemUNIFAC model (dashed line). Our model is shown to be able to 
ft predict the overall composition and the distribution of polymers at high concentrations of formaldehyde. At 383K, 

78 the pressure prediction (c) was not as accurate as composition. Nonetheless, the figure 5 shows that the higher the 

79 equilibrium temperature, the more accurate the calculation of the pressure. Also, chemUNIQUAC model gives a 

80 better description of the azoetropic point. 
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Figure 3: Prediction of formaldehyde-water system vapor-liquid and chemical equilibria at 383K: (a) overall composition of formaldehyde in vapor 
phase vs overall composition of formaldehyde in liquid phase, (b) concentration of methylene glycol (MG) and polyoxymethylene glycols (MG2, 
MG3, MG4) in chemical equilibria and (c) Pxy isothermal phase diagram, (x) experimental values from the literature [17]. Solid line: predicted 
phase diagram with chemUNIQUAC. Dashed line: predicted phase diagram with chemUNIFAC. 



Figure 4: Prediction of isothermal phase diagrams for water-formaldehyde system at four different temperatures (a) 363K (b) 383K (c) 413K and 
(d) 423K. (x) experimental values from the literature [17]. Solid line: predicted phase diagram with chemUNIQUAC. Dashed line: predicted phase 
diagram with chemUNIFAC. 
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Figure 5: Prediction of Pxy isothermal phase diagrams for water-formaldehyde system at four different temperatures (a) 363K (b) 383K (c) 413K 
and (d) 423K. (x) experimental values from the literature [17]. Solid line: predicted phase diagram with chemUNIQUAC. Dashed line: predicted 
phase diagram with chemUNIFAC. 



Figure 6: Distribution of methylene glycol (MG) and polyoxymethylene glycols (MG2, MG3, MG4) in chemical equilibria at different temper¬ 
atures: (a) 338K, (b)353K, (c) 368K and (d) 383K. (x) experimental values from the literature [17]. Solid line: predicted phase diagram with 
chemUNIQUAC. Dashed line: predicted phase diagram with chemUNIFAC. 
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Table 8: Average overall deviation of the gas-phase composition and average overall deviation of the pressure for the water-formaldehyde system 
using chemUNIQUAC model and chemUNIFAC model. 


A y FA (%) A P(%) 

Temperature (K) chemUNIQUAC chemUNIFAC chemUNIQUAC chemUNIFAC 


363 

383 

413 

423 


4.41 4.93 

2.89 3.72 

0.86 1.70 

1.75 3.89 


2.19 2.20 

1.51 1.46 

0.60 0.96 

0.36 1.53 


Table 9: Average deviation in liquid-phase of the polymer distribution for the water-formaldehyde system using chemUNIQUAC model and 
chemUNIFA C model (1/2) _ 

A xmg Ax MGi 

Temperature (K) chemUNIQUAC chemUNIFAC chemUNIQUAC chemUNIFAC 


338 

353 

368 

383 


1.59 6.97 

1.82 6.85 

1.28 7.11 

3.05 8.33 


5.89 13.3 

5.64 12.7 

6.09 13.6 

5.60 11.9 


Table 10: Average deviation in liquid-phase of the polymer distribution for the water-formaldehyde system using chemUNIQUAC model and 
chemUNIFA C model (2/2) _ 


Temperature (K) 

Axmg, 

chemUNIQUAC chemUNIFAC 

A xmGa 

chemUNIQUAC chemUNIFAC 

338 

5.08 

15.2 

5.25 

18.1 

353 

4.62 

13.5 

6.05 

17.7 

368 

6.43 

11.5 

11.5 

14.5 

383 

10.7 

18.0 

19.0 

26.8 


181 Figures 4 and 5 show results for formaldehyde-water system at four temperatures: (a) 363K, (b) 383K, (c) 413K, 

182 (d) 423K. Figure 4 presents the overall composition of formaldehyde in vapor phase vs overall composition of 

183 formaldehyde in liquid phase. Good agreement was obtained between experimental data and the prediction of the 

184 overall composition with the chemUNIQUAC model. Figure 4 confirms that the chemUNIQUAC model provides a 
las better description of the composition at high concentration of formaldehyde. The figure 5 shows the Pxy isothermal 
las phase diagram where the higher the temperature, the better the prediction of pressure. 

is? Figure 6 shows the distribution of methylene glycol (MG) and polyoxymethylene glycols (MG2, MG3, MG4) at 

188 four temperatures: (a) 338K (b) 353K (c) 368K and (d) 383K. This figure points out the interest of the chemUNIQUAC 

189 model which gives good prediction of the polymers distribution in the liquid phase. 

190 4.2.3. Formaldehyde-Methanol System 

191 Table 11 compares the average overall deviation of the gas-phase composition and average deviation of the pres- 

192 sure for the methanol-formaldehyde system between experimental data from the literature [23] and the two models 

193 (chemUNIQUAC and chemUNIFAC). 

194 The values of the deviations reported in table 11 indicate that the proposed chemUNIQUAC model provides a 

195 good representation of the formaldehyde-methanol system, equivalent to the chemUNIFAC. Figures 7 to 9 show the 

196 prediction of formaldehyde-methanol system vapor-liquid and chemical equilibria at three temperatures (333K, 343K, 

197 and 353K). The part (a) of each figure illustrates the overall composition of formaldehyde in vapor phase vs overall 

198 composition of formaldehyde in liquid phase and the part (b) the Pxy isothermal phase diagram. Good agreement with 

199 experimental data was obtained with both models. Experimental distributions of the polyoxymethylene hemiformals 
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Table 11: Average overall deviation in gas-phase of the composition and average overall deviation of the pressure for the methanol-formaldehyde 
system using chemUNIQUAC model and chemUNIFAC model. 


Temperature (K) 

A y FA {%) 

chemUNIQUAC chemUNIFAC 

A P(%) 

chemUNIQUAC chemUNIFAC 

333 

9.08 

6.33 

2.65 

1.89 

343 

5.04 

4.67 

2.57 

3.72 

353 

5.26 

5.16 

3.21 

4.21 


could not be plotted because not available in the literature. Indeed, formaldehyde is rarely stored with pure methanol; 
aqueous solutions of formaldehyde with or without methanol are the usual way to handle it. 

4.2.4. Formaldehyde-Water-Methanol System 

Table 12 compares the average overall deviation of the gas-phase composition and average deviation of the pres¬ 
sure for the formaldehyde-water-methanol system between experimental data from the literature [18] and the two 
models (chemUNIQUAC and chemUNIFAC), calculated with equations 22 to 24. 

The values of the deviations show that the chemUNIQUAC model provides more accurate description of the 
formaldehyde-water-methanol system than the chemUNIFAC model. 


Table 12: Average overall deviation in gas-phase of the composition and average overall deviation of the pressure for themethanol-formaldehyde 
system using chemUNIQUAC model and chemUNIFAC model. 

AP(%) 

Temperature (K) chemUNIQUAC chemUNIFAC chemUNIQUAC chemUNIFAC chemUNIQUAC chemUNIFAC 

333 0.66 0.70 0.93 4.07 0.25 0.26 


Table 13 reports a comparison between experimental data in the literature [18], predicted equilibria with chemUNI¬ 
QUAC and with chemUNIFAC. Overall deviation in the gas-phase and deviation of the pressure between experimental 
data and models were calculated for each equilibrium data. 

In this case also, the values of the deviations for each composition confirm that the chemUNIQUAC model pro¬ 
vides a better representation of the formaldehyde-water-methanol system than the chemUNIFAC model. 

5. Conclusions and Perspectives 

The improved approach using UNIQUAC coupled with chemical reactions (so-called here as the chemUNIQUAC 
model) was proved to quantitatively describes the reactive vapor-liquid equilibrium of formaldehyde-water system, 
formaldehyde-methanol system and formaldehyde-water-methanol system. The chemUNIQUAC model was shown 
to provide a better description of the composition at high concentration of formaldehyde than the original UNIFAC 
model. This model provides a reliable description of the vapor-liquid and chemical equilibria from about 330K 
to 420K, for dilute as well as for concentrated solutions. It also accurately describes the oligomer distribution in 
the liquid-phase. The main interest of our model lies in its good ability to predict the overall composition and the 
distribution of polymers at high concentrations of formaldehyde. 
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Figure 7: Prediction of formaldehyde-methanol system vapor-liquid and chemical equilibria at 333K: (a) overall composition of formaldehyde in 
vapor phase vs overall composition of formaldehyde in liquid phase, (b) Pxy isothermal phase diagram, (x) experimental values from the literature 
[23]. Solid line: predicted phase diagram with chemUNIQUAC. Dashed line: predicted phase diagram with chemUNIFAC. 
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Figure 8: Prediction of formaldehyde-methanol system vapor-liquid and chemical equilibria at 343K: (a) overall composition of formaldehyde in 
vapor phase vs overall composition of formaldehyde in liquid phase, (b) Pxy isothermal phase diagram, (x) experimental values from the literature 
[23]. Solid line: predicted phase diagram with chemUNIQUAC. Dashed line: predicted phase diagram with chemUNIFAC. 




Figure 9: Prediction of formaldehyde-methanol system vapor-liquid and chemical equilibria at 353K: (a) overall composition of formaldehyde in 
vapor phase vs overall composition of formaldehyde in liquid phase, (b) Pxy isothermal phase diagram, (x) experimental values from the literature 
[23]. Solid line: predicted phase diagram with chemUNIQUAC. Dashed line: predicted phase diagram with chemUNIFAC. 
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